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Cli! Abstract 

^, 

'^ , Higher order quantum effects on the magnetic phase diagram induced by four-spin ring exchange 

-(— > , 

on plaquettes are investigated for a two-dimensional quantum antiferromagnet with S = 1/2. 
_L ' Spatial anisotropy and frustration are allowed for. Using a perturbative spin-wave expansion up 

Q , to second order in 1/5 we obtain the spin-wave energy dispersion, sublattice magnetization, and 

CJ ■ 

the magnetic phase diagram. We find that for substantial four-spin ring exchange the quantum 

fluctuations are stronger than in the standard Heisenberg model. A moderate amount of four-spin 
ring exchange couplings stabilizes the ordered antiferromagnetic Neel state while a large amount 



(N 
> 

^Sj , renders it unstable. Comparison with inelastic neutron scattering data points toward a moderate 



ring exchange coupling of 27% to 29% of the nearest-neighbor exchange coupling. 
PACS numbers: 75.10.Jm, 75.40.Mg, 75.50.Ee, 73.43.Nq 



I. INTRODUCTION 

Despite the intense experimental and theoretical activities to understand the origin of 
high temperature superconductivity in layered oxide high-temperature superconductors, the 
underlying microscopic mechanism is still incomplete.-"- Very recently the crucial role of 
magnetic excitations in these compounds has been supported by their observation in the 
whole Brillouin zone up to high energies and high levels of doping.— 

The conventional route to theoretically investigate the magnetic properties of these un- 
doped compounds is the two-dimensional (2D) antiferromagnetic (AF) spin-1/2 Heisenberg 
model with nearest neighbor (NN) AF coupling Ji and next-nearest neighbor (NNN) anti- 
ferromagnetic coupling J2.— For concreteness, we give the studied Heisenberg Hamiltonian 
for a S* = 1/2 antiferromagnet on a square lattice 

i i i 

+ '^^2^ [(Si ■ Si+5^)(Si+5^ ■ Si+s^+Sy) + (Si ■ Si+5j(Si+5^ ■ Si+s^+Sy) 
i 

- (Si ■ Si+5^+5j(Sj+5^ ■ Si+5^) . (1) 

We consider four different exchange couplings: Ji for nearest neighbors (NN) along the rows, 
J[ for NN along the columns, J2 for the next nearest neighbors (NNN) along the diagonals, 
and finally the four-spin ring exchange interaction K. All interactions are assumed to be 
antiferromagnetic, i.e., Ji, J(, J2, K > Q. Here i runs over A^ lattice sites and 5x, Sy are unit 
vectors in both directions. In the present work, we study the parameter region where the 
ground state is of Neel type as shown in Fig. [H We take Ji as the fundamental energy scale 
so that the ground state and its properties depend on the dimensionless ratio 77 := J2/J1 
parametrizing the degree of frustration, the ratio C := J[/Ji parametrizing the degree of 
spatial anisotropy, and the ratio fi = KS'^/Ji parametrizing the relative strength of the 
four-spin ring exchange. Note that the full cyclic permutation around a plaquette comprises 
also two-point couplings along the plaquette edges and along the diagonals.— But they do 
not need to be considered separately because they are incorporated in Ji, J[, and J2- 

Experimentally the ground state phase diagram of these frustrated spin systems can be 
explored from high values to low values of rj by applying high pressures. For example. X-ray 
diffraction measurements on Li2VOSi04 show that the value of rj decreases by about 40% 




AF Phase 

FIG. 1: Classical antiferromagnetic ground state (Neel state) and the various couplings: Ji, J[ 
are nearest neighbor interactions along the row and column directions respectively, J2 is the next- 
nearest neighbor interaction along the diagonals, and K is the cyclic four-spin ring exchange. All 
couplings are assumed to be antiferromagnetic, i.e., Ji, J{, J2, K > 0. 



with increase in pressure from zero to 7.6 GPa.^'^- 

Theoretically, evidence for sizable four-spin ring exchange K in high-temperature 
superconductors^^"— was found soon after the discovery of high-temperature 
superconductivity.— Such exchange processes turned out to be the dominant sublead- 
ing correction to the NN Heisenberg Hamiltonian if it is derived from a three-band 
Hubbard modeU^"— or from a single-band Hubbard model^i^— . Experimental evidence for 
ring exchange stems from the analysis of infrared absorption^^, of Raman response^^i^, 
and of inelastic neutron scatteringi'^'^^i^E. The results indicate that the ring exchange 
coupling reaches between Xring = 2K/Ji = 0.2 and 0.25 relative to the NN coupling. Note 
that for S = 1/2 one has x^ing = 8fi. These findings and the quantitative estimates are 
strongly supported by the analysis of two-leg spin ladder systems such as Sri4Cu2404i, 
Ca8La6Cu2404i, and (Ca, La)i4Cu2404i.— 1^1^"— 

The recent discovery of superconductivity in the class of iron pnictide has ushered a re- 
newed interest in this exciting field.— The parent phases of these materials have been found 
to be metallic, but with columnar AF order.—"—. Since the superconductivity appears in 
immediate proximity of the magnetically ordered phase, it is evident that the magnetic ex- 
citations play an important role.—"— Neglecting the metallicity of the parent phases the 
magnetic excitations can be described by frustrated two-dimensional Heisenberg Hamilton!- 



ans with S > l/2M-~— although the three-dimensionahty cannot be neglected^"—. Ab initio 
calculations seem to indicate a strong spatial anisotropy C ~ of the NN couplings^ fitting 
to the experimental findings.—"— But the weak structural distortion does not explain this 
strong anisotropy. So either orbital order— i^ or higher order magnetic exchange such as NN 
biquadratic coupling^"— may effectively explain the anisotropy. 

Another class of magnetic materials described by the Hamiltonian in Eq. ([1]) are vana- 
dium phosphates. Extensive band structure calculations^ yielded four different exchange 
couplings: Ji and J'^^ between the NN and J2 and Jg between NNN in the compounds 
Pb2VO(P04)2, SrZnVO(P04)2, BaZnVO(P04)2, and BaCdVO(P04)2. For example C ~ 0.7 
and J2/J2 ~ 0.4 were obtained for SrZnVO(P04)2. Also the compound (NO)Cu(N03)3 pos- 
sibly realizes the Ji-J[-J2 model.— 

The above examples corroborate the relevance of the model ([1]). 

It is now well known that at low temperatures the spin- 1/2 antiferromagnetic J1-J2 model 
on a square lattice exhibits new types of magnetic order and novel quantum phases.— For 
J2 = and K = the ground state is Neel ordered at zero temperature. Addition of next- 
nearest neighbor (NNN) interactions induces a strong frustration and breaks the Neel order 
at a quantum critical point J2/J1 ~ 0.4 as found by 1/S expansions^"—, series expansion 
about the Ising limit^, and the coupled-cluster approach^. We stress that the precise nature 
of the phase beyond the Neel phase is still intensely debated^"—. 

A generalization of the frustrated J1-J2 model is the Ji-J[-J2 model where ( = J[/Ji is 
the directional anisotropy parameter.— "^i^ Recently, the role of directional anisotropy on 
the magnetic phase diagram has been investigated in detail using a spin-wave expansion.— 

The next generalization consists in the inclusion of the four-spin ring exchange interaction 
K which is the next important coupling after the NN exchange coupling. Using linear 
spin- wave theory its effects on the magnetic properties of the J1-J2-K model were studied 



in Ref. 



16 



28 



where a quasiclassical phase diagram in 0{{1/S)^) was obtained. In Ref. 
corrections to the spin-wave spectrum to first order in 1/S were studied for finite K using self- 
consistent spin-wave theory. The self-consistent spin-wave theory is a mean-field approach 
which captures only a part of the second-order effects 0{{1/Sy) in the phase diagram. In 
particular, it does not take virtual excitations of two and four magnons into account. To 
consider them a perturbative spin- wave expansion up to l/S"^ is needed. That is the goal of 
the present work. 



In the present paper we investigate the higher-order quantum corrections due to the 
presence of plaquette four-spin ring interactions on the antiferromagnetic phase diagram of 
the Ji-J[-J2-K Heisenberg model on a square lattice, cf. Eq. ([T]). Our calculations use the 
Dyson-Maleev spin representation which facilitates the calculation significantly compared 
to the Holstein-Primakov representation. The concomitant formalism is presented in the 
next section. Results for the spin-wave energies and the magnetizations of the system are 
presented and discussed in Section lllli A quantitative comparison with experimental data 
is also included. Section IIVI contains a brief summary of our results. 

II. FORMALISM 

Quantum fluctuations play a significant role in the magnetic phase diagram of the system 
at zero temperature. We will investigate the role of quantum fluctuations on the stability 
of the Neel phase. We first express the fluctuations around the classical antiferromagnetic 
ground state in terms of the boson operators using the Dyson-Maleev representation. The 
quadratic term in boson operators corresponds to the linear spin-wave theory, whereas the 
higher-order terms represent spin-wave interactions and virtual processes. We keep terms 
up to second order in 1/5. In the next step we calculate the renormalized magnon Green's 
functions and self-energies. Finally, we calculate the magnon energy dispersion and the 
sublattice magnetization up to and including terms of order l/S"^. 

For the Neel ordered phase NN couplings interact between the A and B sublattices while 
NNN couplings link A and A sites or B and B sites, respectively. The Hamiltonian in Eq. ([1]) 
takes the form 

i i i 

+ 2KJ2 [(Sf ■ Sf )(St ■ Sf) + {Sf ■ S^)(Sf ■ Sf ) - (Sf • S^)(Sf ■ Sf )] , (2) 

i 

where j = i + Sx, k = i + Sx + Sy, i = i + Sy as shown in Fig. [H Beside the directional 
anisotropy parameter ( = J'l/Ji, the magnetic frustration between the NN and NNN spins 
rj = J2/J1, and the cyclic four-spin exchange interaction term /i = KS'^/Ji we use z = 2 for 
the coordination number. This spin Hamiltonian is mapped onto an equivalent Hamiltonian 
of interacting bosons by expressing the spin operators in terms of bosonic creation and 



annihilation operators a\a for "up" sites on sublattice A and b'^,b for "down" sites on 
sublattice B using the Dyson-Maleev representation 



^Ai 



S 



Bj 



2S 
2S 



b] 



^2SY 



J/jU-j, 



(3a) 



(25) 



S^^ = V2Sbj, S%^ = -S + b]bj. (3b) 

Substituting Eqs. ([3]) into ([2]) we expand the Hamiltonian perturbatively in powers of 
1/S as 

H = H^i + Ho + Hi + H2 + --- , (4) 

where Hm is of order l/S*"*"^. Note that H_i is just a number representing the classical 
energy. We do not discuss it further because it is irrelevant for the quantum fluctuations. 
Hence the 1/S expansion will be performed around the unperturbed Hamiltonian Hq which 
is the zeroth order Hamiltonian in this sense. Relative to Hq the terms Hi and H2 are first 
and second order terms, respectively. 

Next the real space Hamiltonian is Fourier transformed to momentum space. Then we 
diagonalize the quadratic part Hq by transforming the operators ak and 6k to magnon 
operators a^ and /3k using the usual Bogoliubov (BG) transformations 

«k = ^k^k + "^k/3-k, &-k = mkOk + ^k/3-k- 



The coefficients /k and mk are defined as 



Ck 



2ek 



1/2 



mk 



-sgn(7k 



Ck 



26, 



1/2 



-Xk/k, Xk = sgn(7k) 



Ck 



l + Ck 



1/2 



with 7fcx = cos{kx), '^ky = cos{ky) and 



^k = (1-7^)^^ 
Tik 
Kk ' 
(1 - Afl)-fkx + (C - 4/i)7fcy 



7k 



7ik 

72k 

«k 



1 + C - 8/i 

7fcx7fcj/5 

2(r/ - 2/i) 



1- 



(l-72k). 



(5) 

(6) 

(7a) 
(7b) 

(7c) 
(7d) 
(7e) 



1 + C - »/w 

The function sgn(7k) keeps track of the sign of 7k in the first Brillouin zone (BZ). After 
these transformations, the quadratic part of the Hamiltonian takes the form 

Ho = JiSz{l + C - 8/i) ^ /Ck (ek - 1) + JiSz{l + ( - 8/i) J^ K^e^ («j,«k + /^l/Sk) • (8) 



The first term is the quantum zero-point energy and the second term represents the excita- 
tion energy of the magnons within linear spin- wave theory (LSWT).— 

The part Hi comprises 1/S contribution to the Hamiltonian. We follow the same proce- 
dure as described above. The resulting expression after transforming the bosonic operators 
to magnon operators is 

J,Sz{l + C - 8/.) ^ J^^ i^^t ^^ ^ ^t^^) + B^ (^^t ^t ^ + ^_^«^) - 

k 

■ ^ 6g{1 + 2 - 3 - 4)/i/2/3^4 V^i2;34aia2Q^3a4 + 2V/2;34ai/3-2a3a4 



Hi 



2S 
JiSz{l + C-8fi) 



2SN 



1234 



2^i?;34«I«2/3l3«4 + ^Vi%A(^3Pl,P-2 + ^vSl^pl^a^P-^P^i + 2vSl,l3l,l3l,all3_i 



K7) ^t^tot flt 



.(8) 



(9) flt «t 



V;'2:34«I«2/3-3/3-4 + V;'2;34/3-l/3-2«3«4 + Vit,,(3L,/3U(3^2(3^i 



(9) 



In the above equation momenta ki, k2, ka, k4 are abbreviated as 1, 2, 3, and 4. The first term 
in Eq. (Q is obtained by normal ordering the products of four boson operators with respect 
to creation and annihilation in the magnon operators, i.e., magnon creation operators appear 
always to the left of magnon annihilation operators. The coefficients A^ and -Bk read 



^k 
5k 



Ai- 



KkCk L 

Bi 7ik 



/^k - 7ik 



A, 



ek 



72k 



A, 



ek L 



(1 + 72k) - 7k(7fcx + Iky) 



KkCk 

where the shorthands 
Ai = 

A2 

A, 

Bi 



72k 



+ A3 



E 



ek 



7p7ip 



ilkx + Iky) - 7k(l + 72k) 



(10a) 
(10b) 



2(r/-4/i) /2 



1 + C - 8/i 

4/i /2 

l + C-8/iVA^ 

2(r/ - 2/i) ^2 
1 + C - 8/i 






72p 



+ 72p-7p(7pa;+7py) 



(I)e£[ 



72p - 7p7ip 



(11a) 
(lib) 
(lie) 
(lid) 



are used. 

The second term in Eq. ([9]) represents scattering between spin-waves where the delta 
function (5g(1 + 2 — 3 — 4) ensures that the momentum is conserved within a reciprocal 
lattice vector G. Explicit forms of the vertex factors V^/^4' ' ' ' are given in Appendix [Bl 



The second order term, H2 is composed of six-boson operators and is only present when 
/i 7^ 0. Before the Fourier and BG transformations H2 is of the following form 

8/i5 



Ho 



{2sy- 



y^ {a\ai + h^jhj + aihj + a|6])(a|.afe6j6^ + -a\akakhi + -alh\h\hi 



+ {alaib^jbj + -a\aiaihj + -a\h^jhihj){a\ak + ftjfe^ + a^fe^ + a\h\) 

+ {a\akh^-hj + -a\akakhj + -a|,6]6]6j)(ajai + 6^6^ + 0^,6^ + affej) 

1 1 

- {a\ai + 4afc - a^al - a\ak){h]hjh\hi - -b]b]bjbe - -bjb\b\bi) 

- {a\aia\ak - -a\aiaia\ - -a\a\akak){b]bj + b\bi - bjb\ - b^-bi) 



(12) 



After Fourier and BG transformations to magnon operators Ok? /5k the Hamiltonian in 
normal-ordered form reduces to 



Ho 



AjizS 
l2Sf 



Y^ fcik (aj^aic + /3i^/3k) + Cak (a^/Slk + /S-k^k 



(13) 



The dotted terms contribute only to higher than second order corrections and are thus 
omitted in our calculations. The coefficients Cik and C2k are given in Appendix O 

The quasiparticle energy E^^ for magnon excitations, measured in units of JiSz{l + C ^ 
8/i) up to second order in 1/ S is given as 



Ef = E, 



-A 



sL^i(k,i^k) ^^ 



2Eu 



(14) 



{2S) " (25)2 

Expressions for the magnon Green's functions and self-energies are given in Appendix \M 
The dynamic contributions to the second order self-energies S*^^) ^^-g second order in the 
vertex factors V^^\ These are the contributions which are missed by self-consistent spin- 
wave theory. 

The sublattice magnetization Maf for the A sublattice can be expressed as 

Ml M2 



Maf = S- (altti) = S - AS 



(25) {2Sy 



(15) 



where 

AS = ^y — --, (16a) 



iV^2ek 2' 

k 



^. = ^E^. (16b) 



k 



N Z^\ ^k^ k;^^2 ^ ^ 

k k K 






' /_^ 2/kipiqik+p- 



k,p,q,[k+p-q] [k+p-q],q,p,k 



^^U <WL.^^..^.K 



pq 



(-Ek + -Ep + -Eq + -Ek+p-q)^ 

(16c) 



_l_ ^ ^ k,p,q,[k+p-q] [k+p-q],q,p,k 

El- (Ep + Eq + i5;k+p_q)2 

The zeroth-order term AS* corresponds to the reduction of magnetization within LSWT, M\ 
term corresponds to the first-order \jS correction, and M2 is the second-order correction. 
Again, the parts which are second order in the vertex factors are not captured by self- 
consistent spin-wave theory. 

III. RESULTS 

1 . Spin- Wave Energy 

We obtain the spin-wave energy 2JiS'(l + C " 8n)E^^ for S = 1/2 as a function of 
momenta {kx, ky) for several values of (, rj, and fi by evaluating Eq. JHM in the first BZ. 
For the numerical summation we divide the first BZ in a mesh of N^ points with A'^^ = 48 
and then the contributions from all the points are summed up to evaluate the third term in 
Eq. (TT4|) . In the Dyson-Maleev formalism, no cancellation of divergences occurs so that the 
convergence of the numerical results for A'^^: — !■ C)0 is very good. This is a crucial advantage 
over the use of the Holstein-Primakov representation. We estimate that the results for 
A'^i;, = 48 will not change more than in the third digit if A'^^: is chosen larger. 

Figure |2] shows a comparison between the results from LSWT (long-dashed lines), first- 
order (dot-dashed lines) and second-order corrections (solid lines) to the spin-wave energy 
spectrum for isotropic coupling (C = 1) for two choices of frustration and and ring exchange. 
For the moderate value /i = 0.025 corresponding to 2K/Ji = 0.2 the 1/S correction is 
substantial while the l/S"^ correction is fairly small. This is very similar to the corrections 
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FIG. 2: (Color online) Spin-wave energy EQ^ / Ji obtained from LSWT (long-dashed lines), with 
1/5 (dot-dashed lines) and with 1/ S"^ corrections (solid lines) for the Neel-ordered phase. We have 
chosen spatially isotropic coupling (" = 1. In the left panel we show the corrections for relative 
frustration 77 = and ring exchange ^ = 0.025; in the right panel we show them for frustration 
T] = 0.2 and ring exchange /i = 0.12. In the latter case, the 1/5^ terms in the Hamiltonian provide 
significant corrections to both the LSWT and 1/5 results. 

for the NN Heisenberg model at // = q.— "— 1^"— The right panel of Fig. [2] tells quite a 
different story. For substantial ring exchange the quantum corrections are very large and 
cannot be ignored. We point out that this is not due to the frustration alone as can be seen 
by inspecting the results for substantial values of r/, but without ring exchange /x = 0, in 



Ref. 



60 



|. The l/S*^ corrections for /i = are as small as they are for the NN Heisenberg 
model, in contrast to the result in the right panel of Fig. [2l 

In the panels of Fig. [3] the evolution of the spin-wave energy spectrum including cor- 
rections up to second-order for various values of (,^7] and /i are shown. The spin-wave 
dispersions for the couplings C, = 1 and C = 0.4 at /i = were reported earlier using 
the Holstein-Primakov representation.— The results from the Dyson-Maleev and from the 
Holstein-Primakov representation coincide as it has to be for physically observable results 
of a systematic expansion in a small parameter. 

For /i = and 77 = 0, the energy at (7r/2,7r/2) is larger than the energy at (7r,0), cf. 
upper left panel in Fig. [3l This dip of the dispersion at (vr, 0) has been first computed by 
high-order series expansion (HSE) around the Ising limit^i^ and was confirmed by quantum 
Monte Carlo calculation (QMC).— HSE and QMC find that the dip is about 9% deep, i.e.. 
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FIG. 3: (Color online) The effect of ^ on the spin-wave energy E^^ / Ji for the Neel-ordered phase 
with 1/5^ corrections is shown for two values of r/ = 0, 0.2 and C = 1- 

[^((7r/2,7r/2)) - E((7r, 0))]/^((7r/2, 7r/2)) ^ 0.09. Experimentally, the dip is found to be 
about 7% in compounds in which no couplings beyond Ji are thought to play a role, in 
reasonable agreement with HSE and QMC-"^ 

In contrast, LSWT and order 1/ S do not find a dip at all. In order l/S*^, it is present 
but as small as 1.4% and in order 1/S''^ it takes the value of 3.2%.— Thus one must be aware 
that the data in Fig. [3] does not capture all aspects of the dispersion between (tt, 0) and 
(7r/2, 7r/2). But in the remaining BZ the significance of corrections of third order and higher 
is rather small and the agreement with the series expansion results very good. 

Having the above minor caveat in mind, we discuss the much stronger influence of frus- 
tration and of ring exchange in the following. Increasing the value of /i to positive values 
the energy at (7r/2,7r/2) decreases more strongly than the one at (vr, 0), see left panels of 
Fig. [3l Hence, beyond some finite value of four-spin ring exchange there is a dip from (vr, 0) 
to (7r/2,7r/2). This agrees qualitatively with experimental findingsii^, which see a 13% dip, 
and with an analysis based on self-consistent spin-wave theory.— Even larger values of /i will 
lead to a complete softening of the magnon mode at (7r/2, vr/2). This indicates a competition 
between an ordered orthogonal state at modulation (7r/2, 7r/2) and the ordered Neel state 
at (tt, tt) upon increasing /x. 
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Another important issue is the effect of finite frustration 77 > which has been investigated 
before without ring exchange.— 1^ Indeed, finite frustration induces a significant dip at (tt, 0) 
relative to (7r/2,7r/2), i.e., £'((7r,0)) < E{{tx /2,tx /2)), so that frustration pushes the system 
into the opposite direction as does the ring exchange. But in the presence of substantial 
ring exchange the effect is reversed: Comparing the upper and lower left panels in Fig. [3] 
and inspecting Fig. |4] we see that increasing frustration supports the tendency to soften the 
mode at (vr/2, vr/2) which will eventually destabilize the Neel order. 

Spatial anisotropy, see right panel in Fig. |3l does not alter this picture qualitatively. A 
strong anisotropy C, <1 seems to support the tendency to mode softening and the concomi- 
tant destabilization of the Neel order. 

2.5 1 



C=i.o, \i=o.n 







(71, 0) (71/2,71/2) (3t,3t) 



(71,0) 



FIG. 4: (Color online) Spin-wave energy E^^ / Ji including 1/5^ corrections for (^ = 1, ^u = 0.12 for 
various values of r/. 



2. Quantitative Analysis of the Inelastic Neutron Scattering Data 

We use our model to quantitatively analyse the experimental data obtained in Ref. 
by inelastic neutron scattering for La2Cu04. We disregard any spatial anisotropy because 
La2Cu04 is tetragonal so that we set C = 1. The experimental data displays a significant 
dip at (7r/2, 7r/2) relative to the energy at (tt, 0). This points toward a sizable four-spin ring 
exchange^i^. 
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FIG. 5: (Color online) Comparison of the measured spin- wave energy E^^ as obtained by inelas- 
tic neutron scattering in La2Cu04 with the theoretical results including 1/5^ corrections for the 
spatially isotropic model {(^ = 1) for Nl = 24. For given moderate values i] of relative frustra- 
tion a value fi of the four-spin ring exchange can be found such that the dispersions match the 
experimental data. 

Our findings are shown in Fig. [51 They strikingly confirm that substantial values of n 
are needed to explain the observed energy dip at (7r/2,7r/2). For instance, for 77 = one 
needs /i = 0.034, and Ji = 143 meV; for 77 = 0.01 /i = 0.036, and Ji = 146 meV; for 
T] = 0.02 fi = 0.0375, and Ji = 148 meV (not shown). Even for 77 = 0.10 the parameters 
/i = 0.046, and Ji = 174 meV yield a theoretical dispersion which is indistinguishable from 
those displayed in Fig. [51 Note that the agreement of the steeply rising parts of the dispersion 
is not completely perfect because the theoretical curves remain a bit below the experimental 
data points. 

We conclude that from the experimental data for the spin-wave energies the relative 
frustration and the relative ring exchange cannot both be determined independently. Based 
on the results of systematic derivations of extended Heisenberg models for the cuprates 
starting from microscopic Hubbard models^^i^^i^i we stick to small values of frustration rj ^ 
0.01. According to our fits this implies Xrmg = 2K/Ji = 8fi = 0.29. This relative four-spin 
ring exchange is slightly larger than we would expect from the systematic derivations.— >^>^ 
It is also slightly larger than the value 0.24 found in the analysis by self-consistent spin-wave 
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theory.— 

On the one hand, the agreement is good in view of the remaining uncertainty in the 
description of the spin-wave energies at the zone boundary, see our discussion above. On 
the other hand, a further improved theoretical treatment of spin-waves is desirable. 

3. Sublattice Magnetization and the Phase Diagram 

We calculate the sublattice magnetization Map from Eq. (1151) by numerically evaluating 
Eqs. fll6al) - fll6cp with C = 1 ^^^ 0-8 ^"^^ fo^^ /^ = 0,0.12 and 0.22. Especially to obtain the 
second order correction term M2 we sum up the values of N^/A points of k in a quarter of 
the first BZ and N^ points of p and q in the first BZ, with A'^^^ = 36 sites along one axis. 

Figure [6] shows the sublattice magnetization with increase in the frustration parameter 
7] = J2/J1 for the isotropic case ( = J[/Ji = 1 for three different values of plaquette ring 
exchange coupling /i = KS^/Ji = 0,0.12, and 0.22. For each case, three different curves 
are plotted: The long-dashed lines represent the LSWT prediction, the dotted lines include 
the first-order (l/S) correction to the LSWT results, and the solid lines include corrections 
up to second-order (1/5^). Upon increasing frustration the dotted curves of the first-order 
corrections diverge. However, 1/S^ corrections (M2) significantly increase with frustration 
and stabilize the apparent divergence of the magnetization. We find that the magnetization 
with second-order corrections decreases steadily at first with increase in rj and then sharply 
drops to zero at a critical value of 77 = rjc- Assuming that the Neel phase loses its stability 
continuously, rjc marks the quantum critical point at which the AF order is destroyed and the 
system enters into another state characterized by other types of order. The precise order 
of the phase transition and the nature of the subsequent phase is still matter of intense 
debate.-^--^. 

Without four-spin ring exchange, i.e., /i = 0, Map with second-order corrections begins 
from 0.307 at 77 = and decreases upon rising frustration till 77 ~ 0.32. Finally it vanishes 



60 



usmg 



at ?7ci ~ 0.411. For this case, we reproduce the magnetization plot obtained in Ref. 
a similar perturbative l/S expansion based on the Holstein-Primakov representation. The 
LSWT prediction for the critical point is lower at ~ 0.38. With increase in the four-spin 
ring exchange /i the values of the magnetization at r/ = increase. For example, we find 
MAp(r/ = 0,/i = 0.12) ^ 0.458 and Maf{v = 0,/i = 0.22) ^ 0.524. These numbers are 
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significantly larger than the predictions from LSWT which are 0.381 and 0.466, respectively. 
We conclude that without NNN frustration (77 = 0) the pure four-spin coupling /i favors 
the Neel order. This is in qualitative accord with the observation that the spin gap of the 
disordered paramagnetic phase of spin ladders is reduced on increasing four-spin coupling 
^.i^i^i^i^ Thus finite four-spin coupling pushes spin ladders closer to a gapless phase which 
is likely to display quasi-long range order with powerlaw correlations. 

We observe that first and second order corrections provide significant contributions to 
the entire magnetization curves. For small //, the corrections M2 start from a small positive 
value and then switch sign and become negative with increase in 77. However, for large n, 
say ji = 0.22 M2, corrections are negative throughout. 
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FIG. 6: (Color online) The sublattice magnetization Maf is plotted for C = 1 and for three different 
values of yLi = (black), 0.12 (blue/dark gray), 0.22 (orange/light gray) as a function of the relative 
magnetic frustration 7]. For all three cases, results from linear spin- wave theory (dashed lines), 
with 1/S (dot-dashed lines), and with 1/S'^ corrections (solid lines) are shown. Magnetization 
curves with 1/S corrections alone diverge in all cases. However, 1/S'^ corrections compensate the 
divergence and the magnetization curves steadily decrease to zero at critical values r]c We find 
rjc = 0.411 (^f = 0), 0.423 (^ = 0.12), and 0.399 (^i = 0.22). 

Another interesting feature portrayed in Fig. [6] is the change in the critical value of rj 
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with fi. For /i = the magnetization vanishes at the critical value of frustration r]c ~ 0.411. 
With increase in /i, the value of rjc increases initially till a turning value oi fi = fit ^ 0.12 
is reached beyond which rjc decreases again. For example, rjc ~ 0.423 for /i = 0.12, but 
rjc ~ 0.399 for /i = 0.22. This implies that the four-spin ring exchange interaction favors the 
Neel order and thus extends the AF region only for small values. Beyond the turning value 
yU. = yUt is reached the ring exchange coupling destabilizes the Neel phase. This is shown in 
the rjc-fi phase diagram in Fig. [3 
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FIG. 7: (Color online) rjc-fj. phase diagram for ^ = 1. With increase in /i, rjc increases up to a 
maximum value 0.423 at /U = ^t ~ 0.12 and then sharply decreases. This shows that the ring 
exchange coupling fj, initially favors the Neel ordering of the NN spins till the turning value fit is 
reached. For fi > //f, the four-spin coupling enhances destabilizes the Neel order. 

Next we study the influence of directional anisotropy between the horizontal and vertical 
NN couplings implying C < 1- This spatial anisotropy does not lead to frustration, but it 
weakens the NN coupling because the vertical NN coupling is lowered. Hence we expect a 
qualitatively similar behavior as before, but at lower values of t] and /i. This expectation is 
confirmed by the following results. 

Figure [H] shows the magnetization upon increasing rj for the spatially anisotropic case. We 
choose ( = 0.4 with the three values of ring exchange coupling /i = 0, 0.08, and 0.13. Here 
the values of the magnetization without NNN frustration are Maf(^ = 0,/x = 0.08) ~ 0.40 
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and Maf{v = 0,/i = 0.13) ~ 0.438. Again these numbers are again larger than the LSWT 
values which are 0.350 and 0.406, respectively. 

It is interesting to observe that with increase in rj the magnetization with just 1/S cor- 
rections (dotted curves) diverge except for the case when /i = 0.13. We find that this 
divergence ceases to occur for /i ^ 0.10. As before, l/S"^ corrections significantly modify the 
magnetization curves. The critical values of t] at which the Neel phase is unstable are 0.176, 
0.191, and 0.15 for yU = 0, 0.08 and 0.13, respectively. The LSWT predictions for these three 
cases are 0.172, 0.188, and 0.194, respectively. Notice that the LSWT prediction rjc = 0.194 
for /i = 0.13 is larger than the value rjc = 0.15 obtained including first and second order 
corrections. 
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FIG. 8: (Color online) Sublattice magnetization Maf with spatial anisotropy C = 0.4 between 
the vertical and the horizontal NN couplings for three values of ^ = (black), 0.08 (blue/dark 
gray), 0.13 (orange/light gray) as a function of frustration r]. For all three cases, results from 
LSWT (dashed lines), with 1/5" (dot-dashed lines), and with l/S*^ corrections (solid lines) are 
shown. Maf with 1/S corrections alone diverge for /i = and 0.08, but not for fj, = 0.13 where it 
converges, cf. main text 



It is worth exploring the infiuence of the spatial anisotropy ( on the T]c-fi phase diagram. 
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This is done in the panels of Fig. |9]for ( = 0.4 and 0.2. The results are qualitatively similar 
to those for (^ = 1 in Fig. [TJ but at lower values of rj and fi as we expected. The Neel phase is 
stabilized by small values of /i. But beyond the turning values Ht the four-spin ring exchange 
starts to reduce the parameter region of the Neel phase. 
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FIG. 9: (Color online) rjc-fJ, phase diagram for ^ = 0.4 (left panel) and 0.2 (right panel), to be 
compared with the phase diagram for the spatially isotropic case C = 1 in Fig. [7] 



IV. CONCLUSIONS 

For S = 1/2 Heisenberg the four-spin ring exchange coupling on plaquettes is the next 
important interaction after the nearest-neighbor exchange. In this work we have investigated 
its influence on the zero temperature magnetic phase diagram of a spatially anisotropic and 
frustrated Heisenberg antiferromagnet on the square lattice. 

In particular, we studied higher-order quantum effects in a systematic perturbative spin- 
wave expansion in the inverse spin 5*. We have calculated the spin-wave energy and the 
magnetization up to and including the second-order corrections. They contribute signifi- 
cantly to the shape of the magnetic phase diagram, especially as the frustration between the 
next-nearest neighbor spins increases. The obtained magnetic phase diagram shows that the 
four-spin ring exchange coupling initially favors the Neel order until a specific turning value 
is reached. Beyond this values a further increase in the ring exchange coupling increases 
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the frustration in the system and reduces the parameter region in which the Neel order 
represents the stable ground state. 

Moreover, we analyzed the available neutron scattering data and found that a ring ex- 
change coupling 2K of about 27% to 29% of the nearest-neighbor exchange is required 
to explain the data. The additional determination of the relative frustration in a three- 
parameter fit is not possible because the dispersions for various triples of nearest-neighbor 
exchange, frustration, and four-spin ring exchange are indistinguishable if the energies at 
(tt, 0) and (7r/2, 7r/2) are matched. 
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Appendix A: Green's functions and Self-energies 

The time-ordered magnon Green's functions are defined as 

G«„(k,t) = -z(TK(t)4(0))), G;3;3(k,t) = -z(T(/3lk(t)/3_k(0))), (Ala) 

G,^(k,t) = -z(TK(t)/3-k(0))), G';3,(k,t) = -z(T(/3lk(t)4(0))), (Alb) 

Considering Hq as the unperturbed Hamiltonian the Fourier transformed unperturbed prop- 
agators are 

G'L(k,a;) = ^— -, G°,(k,a;) = L^, (A2a) 

G%{k,co) = G°„(k,a;) = 0, (A2b) 

with 6 — !■ 0+. The spin- wave energy E'k = Kk^k is measured in units oi JiSz{l + (^ — 8/i). 
The graphical representations of the Green functions are shown in Fig. [TOT a). Note the 
differing convention for the arrows which help to represent the conservation of the total S^ 
component in the diagrams efficiently, see Fig. ^U\ 
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FIG. 10: (a) The solid and the dashed hnes correspond to the a and /? propagators. Second-order 

(2) (2) 

diagrams for the self-energies SQ,a(k, w) and S^i(k, w) are shown in (b) and (c). The diagrams 
in (d) contribute only to T,aaO^,^)- V^"^' ,V^^\V^^' jV^"^' ,V^^' are the vertex factors, see main 
text. Note that at each vertex two arrows enter the vertex and two leave it which reflects the 
conservation of the total Sz component. 



The full propagators Gij{\i.,uj) satisfy the matrix Dyson equation 



UJ 



where the self-energy Sjj(k) can be expressed in powers of 1/(25) as 
S,,(k,a;) = J^S«(k,a;) + ^sg)(k,a;) + . . . . 



(25)- 
The first-order self-energy terms read 

-'afS 



{2sy 






S2(k,c^) = S«(k,a;) = i?, 



(A3) 



(A4) 

(A5a) 
(A5b) 
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The second-order self-energy terms originate from the Feynman diagrams in Figs. [TOT b). 
(c), and (d). The coefficients Cik and C2k stem from the normal-ordering oi'H2- The complete 
expressions read 

P 

K2) 



Si^i(k,cu) 



+ 



pq 

y{r) y(S) 

k,p,q,[k+p-q] [k+p-q],q,p,k 



P 

(3) 



k,p,q,[k+p-q] [k+p-q],q,p,k 



uj — Er, — E„ — E\ 



[k+p-q] 



+ 15 



^(2), 



k,cu) 



w + ^p + 



-Eq + -^[k+p-q] - "^^ 

pq 



C 



2k 



2(^pVkAk+p-q])' X 



V, 



(2) 



v; 



(7) 



k,p,q, [k+p-q] [k+p-q] ,q,p,k 



T/(7) T.(5) 

k,p,q,[k+p-q] [k+p-q],q,p,k 



W — -Ep — -Eq 



i5. 



(A6a) 



(A6b) 



-E'[k+p-q] + ^^ W + -Ep + -Eq + -Efk+p-q] 

where [k + p — q] is meant to be mapped to (k + p — q) in the first BZ by an appropriate 
reciprocal vector G. In deriving Eqs. flA6aP and flA6bp we have used the symmetry properties 
of the vertices, see Eq. (lB2p . 
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Appendix B: Vertex factors 

The expressions for the vertex factors are very lengthy. It is convenient to first define the 
following functions 



Jx 


= 72(1 - 4) + 72(2 - 4) + 72(1 - 3) + 72(2 - 3) - 


72(1) -72(2) -72(1 -3-4) 




- 72(2-3-4), 






J1 


= 72(1 - 4) + 72(2 - 4) + 72(1 - 3) + 72(2 - 3), 




Si 


= 7.(4)7,(2 - 4) + 7,.(1 + 2 - 4)7,(1 - 


-4) + 7.(l- 


-3)7,(1 + 2 -3) +7.(2 -3)7,(3) 




+ 7x(3)7,(2-3) + 7.(l + 2-3)7,(l- 


-3) + 7.(l- 


-4)7,(1 + 2-4) +7.(2-4)7,(4), 


S2 


= 7,(4)7,(1 -4) + 7.(1 + 2 -4)7,(2- 


-4) + 7.(2- 


-3)7,(1 + 2 -3) +7.(1 -3)7,(3) 




+ 7.(3)7,(1 - 3) + 7.(1 + 2 - 3)7,(2 - 


-3) + 7.(2- 


-4)7,(1 + 2 -4) +7.(1 -4)7,(4), 


Ss 


= 7.(1 -3 -4)7,(2 -4) +7.(1)7,(2- 


-3)+7.(l- 


-4)7,(2 -3 -4) +7.(1 -3)7,(2) 




+ 7.(2 -3-4)7,(1 -4) +7.(2)7,(1- 


-3)+7.(2- 


-4)7,(1 -3 -4) +7.(2 -3)7,(1), 


s. 


= 7.(1 - 3 - 4)7,(2 - 3) + 7.(1)7,(2 - 


-4)+7.(l- 


-3)7,(2 -3 -4) +7.(1 -4)7,(2) 




+ 7.(2 -3-4)7,(1 -3) +7.(2)7,(1- 


-4)+7.(2- 


-3)7,(1 -3 -4) +7.(2 -4)7,(1), 


s. 


= 7,(2)7,(2 -3) + 7.(2 -3 -4)7,(2- 


-4)+7.(l- 


-4)7,(1 -3 -4) +7.(1 -3)7,(1) 




+ 7,(1)7^(1 -3) + 7.(1 -3 -4)7,(1- 


-4)+7.(2- 


-4)7,(2 -3 -4) +7.(2 -3)7,(2), 


Se 


= 7,(2)7,(2 -4) + 7.(2 -3 -4)7,(2- 


-3)+7.(l- 


-3)7,(1 -3 -4) +7.(1 -4)7,(1) 




+ 7,(1)7^(1 -4) + 7,.(1- 3 -4)7,(1- 


-3)+7.(2- 


-3)7,(2 -3 -4) +7.(2 -4)7,(2), 


S7 


= 7.(1 + 2 - 3)7,(1 - 4) + 7.(3)7,(2 - 


-4) + 7.(l- 


-3)7,(1 + 2-4) +7.(2-3)7,(4) 




+ 7,(1 + 2 - 4)7,(1 - 3) + 7.(4)7,(2 - 


-3) + 7.(l- 


-4)7,(1 + 2-3) +7.(2-4)7,(3), 


Ss 


= 7.(1 + 2 - 3)7,(2 - 4) + 7.(3)7,(1 - 


-4) + 7.(2- 


-3)7,(1 + 2 -4) +7.(1 -3)7,(4) 




+ 7,.(1 + 2 - 4)7,(2 - 3) + 7.(4)7,(1 - 


-3) + 7.(2- 


-4)7,(1 + 2 -3) +7.(1 -4)7,(3), 


S9 


= 7.(1 -4)7,(2 -4) +7.(1 -3)7,(2 - 


-3)+7.(2- 


-4)7,(1 -4) +7.(2 -3)7,(1 -3), 


<Sio 


= 7,(2 -3)7,(2 -4) +7.(1 -3)7,(1 - 


- 4) +7.(2- 


-4)7,(2 -3) +7.(1 -4)7,(1 -3), 


Sn 


= 7.(2)7,(4) + 7.(4)7,(2) + 7.(1 + 2 - 


- 3)7,(1 - 3 


-4)+7.(l-3-4)7,(l + 2-3) 




+ 7,(1 + 2 -4)7,(1) +7,(1)7,(1 + 2- 


-4)+7.(2- 


-3 -4)7,(3) +7.(3)7,(2-3-4), 


Su 


= 7.(2)7,(3) + 7.(3)7,(2) + 7.(1 + 2 - 


-4)7,(1-3 


-4)+7,(l-3-4)7,(l + 2-4) 




+ 7.(1 + 2 -3)7,(1) +7.(1)7,(1 + 2 - 


-3)+7.(2- 


-3 -4)7,(4) +7.(4)7,(2 -3 -4), 
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7,(1)7,(4) + 7x(4)72;(l) + 7x(l + 2 - 3)7,(2 - 3 - 4) + 7,(2 - 3 - 4)7,(1 + 2-3) 



+ 7,(1 + 2 - 4)7,(2) + 7.(2)7,(1 + 2 - 4) + 7,(1 - 3 - 4)7,(3) + 7.(3)7,(1 - 3 - 4), 
5i4 = 7.(1)7,(3) + 7.(3)7,(1) + 7.(1 + 2 - 4)7,(2 - 3 - 4) + 7,(2 - 3 - 4)7,(1 + 2-4) 
+ 7.(1 + 2 - 3)7,(2) + 7.(2)7,(1 + 2 - 3) + 7.(1 - 3 - 4)7,(4) + 7.(4)7,(1 - 3 - 4). 

The vertex factors required for our calculations are 



V, 



(2) 



12:34 



+ 



+ 



- X37i(2 - 3) - X47i(2 - 4) - xiX2X37i(l - 3) - xiX2X47i(l - 4) 

a;ia;27i(l) + 7i(2) + xiX2X3X47i(l - 3 - 4) + X3X47i(2 - 3 - 4) 
r/-2/i 



1 + C - 8/i 
4/i 



X2 + ^gXiXsXa 



Jl 



(X2 + ^0X3X4,) J2 + -(iSi + X1X2S2 + 0:1X353 + X1X4S4, 



l + c-812 

X2X3O5 + X2X4CIQ + X3X4O7 + XiX23;3X408 — 2Xii_>g — 2X2X3X4010 



X41S11 



V, 



(3) 



12;34 



+ 



3:3512 — X1X2X4O13 — X1X2X3O14J 

- xi7i(l - 3) - X27i(2 - 3) - xiX3X47i(l - 4) - X2X3X47i(2 - 4) 

XiX37i(l) + X2X37i(2) + xiX47i(l - 3 - 4) + X2X47i(2 - 3 - 4) 
r]-2fi 



(Bla) 



1 + C - 8;U 
4/i 



X3 + $gXiX2X4 



Jl 



(X3 + ^gXiX2X4)J2 + -(x2X35i + X1X31S2 + X1X21S3 + X1X2X3X454 



V, 



(5) 



12;34 



+ 
+ 

+ 



I + C-812 
1S5 + X3X41S6 + X2X41S7 + X1X41S8 — 2X1X2X3^9 — 2X41S10 — X2X3X4511 

X21S12 — X1X3X41S13 — X11S14) 

- X2X3X47i(l - 3) - xiX3X47i(2 - 3) - xi7i(2 - 4) - X27i(l - 4) 
xiX47i(2) + X2X47i(l) + xiX37i(2 - 3 - 4) + X2X37i(l - 3 - 4) 



(Bib) 



^-2/x 

1 + C - 8/i 
4/i 



X1X2X4 + ^0X3 



Jl 



(X1X2X4 + $^3^3)^2 + ^(a:iX4iSi + X2X41S2 + X3X41S3 + ^4 



1 + C - 8/i 
X1X2X3X41S5 + X1X21S6 + X1X31S7 + X2X31S8 — 2X41S9 — 2X1X2X31S10 



xiiSii 



- X1X3X4512 - X25i3 - X2X3X4514) 



(Blc) 
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T/(7) 
*^12;34 



a;ia;47i(l - 3) + xiX37i(l - 4) + X2a:37i(2 - 4) + 0:2X471 (2 - 3) 
xiX3X47i(l) - 3:2X3X471(2) - xi7i(l - 3 - 4) - X27i(2 - 3 - 4) 



+ 



1 + C - 8/i 
4/i 



X3X4 - ^gXiX2 



Jl 



(X3X4 + ^0X1X2) J2 + ^(-X2X3X4iSi - X1X3X452 - X1X2X41S3 



1 + C - 8/i 
X1X2X354 - X455 - X356 - X21S7 - Xi^s + 2X1X2X3X4^9 + 2Sio + X2X3S11 



+ X2X4S12 + XiX35i3 + XiX45: 
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(Bid) 



K 



(8) 



12;34 



xiX47i(2 - 4) + X2X47i(l - 4) + xiX37i(2 - 3) + X2X37i(l - 3) 
Xi7i(2) - X27i(l) - XiX3X47i(2 - 3 - 4) - X2X3X47i(l - 3 - 4) 

7/-2/i 



1 + C - 8/i 
4/i 



X1X2 - ^0X3X4 



Jl 



X1X2 + <I>gX3X4) ^^2 + ^(-^I'^l - ^2^52 - X353 - X454 



.l + C-8/i 

— X1X2X355 — X1X2X456 — X1X3X457 — X2X3X45J 
+ XiX45ii + XiX35i2 + X2X4513 + X2X3514) 



2^9 + 2X1X2X3X45 



10 



(Ble) 



where $g = exp(iG^.), Gx being the x-component of the reciprocal lattice vector G appearing 
in the momentum conserving delta- function in Eq. ([9]). These vertex factors fulfill the 
following symmetry relations 



V, 



(2) 



V, 



12;34 

(7) 



K2) 



K3) 



12;34 



^12;43i ''l2;34 
'^21;34 ~ ''12;43j 



yi3) . 
''21;34i 



V, 



(8) 



12;34 



T/(5) _T^(5) 
*^12;34 ~ ''21;345 

T/{8) _y(S) 
*^21;34 ~ *^12;43- 



(B2a) 
(B2b) 



If no reciprocal lattice vector is involved in the momentum conservation, i.e., G = 0, there 
are some additional symmetries 



y(3) _yi5) . 
^12:34 ~ ''l2;34) 



'■^12;34 ~ ''l2;34- 



(B2c) 
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Appendix C: Coefficients Cik and C2k 



We define the functions Pk and Qk 
2 \2 



Pk = (^) J22ilil\xlxl{Q + Q^2{k) + 6^2{2) + 2^2ik-2)+rA^-2)^y{l + 2) 



+ 
+ 
+ 
+ 



lx{k 
lx{k 



12 
1 

1 



2)7,(A; - 1 + 2) + a;f 272(A; - 2) + 672(2) + i^k + 1 - 2)7,(A; 



2) 



2)7,y(/e + 1 - 2) + 7,(A; - 1 + 2)7^(A; 



2)+7,(A;-l + 2)7j,(A;+l-2) 



+ 



Qv 



xiX2[A^x{k - lhy{k - 2) + A-f^k - 2)^y{k - 1) + 47,(l)7j,(2) + 47,(2)7^,(1) 
67,(1 - 2) + 67,(1 - 2) + 47,(A;)7j,(A; - 1 - 2) + 47,(A: - 1 - 2)7,(A;)} 
xlx2[^l.{k)ly{k - 2) + 87,(A: - 2)-iy{k) + 127,(2) + 127^,(2) + 47,(1)7,(1 - 2) 
47,(1 - 2)7^,(1) + 27,(A; - 1)7,(A; - 1 - 2) + 27,(fc - 1 - 2)^y{k - 1) 
27,(A; + 1 - 2)7,(A; - 1) + 27,(A; - 1)7,(A; + 1-2)} 
xi|47,(2)7,(l - 2) + 47,(1 - 2)7,(2) + 27,(A; - 2)7,(A; -1 + 2) 
27,(A; - 1 + 2)7,(A; - 2) + 27,(A; + 1 - 2)7,(A; - 2) + 2y,{k - 2)-fy{k + 1-2) 
7,(1 - 2)7,(1 + 2) + 7,(A; - 1 - 2)7,(A; + 1-2)} . 

J2 2^?^2 [xlxl^G-f^k) + 6-fy{k) + 47,(2)7,(A; - 2) + A^,{k - 2)7,(2) 
12 

7,(fc - 1 - 2)7,(1 - 2) + 7,(1 - 2)7,(A; - 1 - 2) 



(Cla) 



nJ 



+ 

+ xi\ 47, (2)7, (A; - 2) + 47, (/c - 2)7,(2) + 7, (A; 



2)7,(1 - 2) 



+ 
+ 

+ 



7,(1 - 2)7,(A: + 1 - 2) + 7,(A; - 1 + 2)7,(1 - 2) + 7,(1 - 2)7,(A; -1 + 2) 
XiX2{87,(2)7,(A; - 1) + 8-f,{k - 1)7,(2) + 67,(A; - 1 + 2) + 67,(A; -1 + 2) 
47,(A;)7,(1 - 2) + 47,(1 - 2)7,(fc)} 

x?a;2{87,(A;)7,(2) + 87,(2)7,(A;) + 127,(A; - 2) + 127,(A; - 2) + 47,(A; - 1)7,(1 
47,(1 - 2)7,(A; - 1) + 47,(A; - 1 - 2)7,(1) + 47,(1)7,(A; - 1 - 2)} 



2) 



- xi 47,(A; - 2)7,(1 - 2) + 47,(1 - 2)7,(A; - 2) + 47,(2)7,(A; 



+ 47, (fc 



2) 



2)7,(2) + 7,(1 - 2)7,(A; - 1 - 2) + 7,(A; - 1 - 2)7,(1 - 2) 



(Clb) 
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Then, the static second-order corrections are given by 



Cik = (4 + m2)Qk + 24mkPk, (C2a) 

C2k = (4 + "^k)^k + 24mkQk. (C2b) 
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